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Abstract. We model in simple terms the angular momentum (J) problem of galaxy 
formation in CDM, and identify the key elements of a scenario that can solve it. The 
buildup of J is modeled via dynamical friction and tidal stripping in mergers. This 
reveals how over-cooling in incoming halos leads to transfer of J from baryons to 
dark matter (DM), in conflict with observations. By incorporating a simple recipe of 
supernova feedback, we match the observed J distribution in disks. Gas removal from 
small incoming halos, which make the low-J component of the product, eliminates 
the low-J baryons. Partial heating and pufnng-up of the gas in larger incoming halos, 
combined with tidal stripping, reduces the J loss of baryons. This implies a higher 
baryonic spin for lower mass halos. The observed low baryonic fraction in dwarf galaxies 
is used to calibrate the characteristic velocity associated with supernova feedback, 
yielding Vfb ~ 100 km s -1 , within the range of theoretical expectations. The model 
then reproduces the observed distribution of spin parameter among dwarf and bright 
galaxies, as well as the J distribution inside these galaxies. This suggests that the 
model captures the main features of a full scenario for resolving the spin crisis. 

1 Introduction 

The 'standard' model of cosmology, CDM, which assumes hierarchical buildup of 
structure, is facing difficulties in explaining observed properties of galaxies, such 
as the number density of dwarfs and the inner density profile of halos. Standing 
out is the angular- momentum problem, that is the apparent failure of the theory, 
via simulations, to reproduce the large sizes of disk galaxies and their structure. 
We make progress by first reproducing the problem via a simple model in which 
the important physical elements are spelled out, and then incorporating in this 
model the key process which may cure the problem — feedback. 

The sizes of disks are commonly linked to the spins of their parent halos as 
measured in N-body simulations [10]. The assumptions that the baryons and 
DM share the same distribution of specific angular momentum j and that the 
baryons conserve their j while contracting to a disk lead to disk sizes comparable 
to those observed. However, simulations that incorporate gas find that most of 
the baryonic j is transfered to the DM, resulting in disk sizes smaller by an order 
of magnitude [e.g. 14,15], and thus leading to a spin catastrophe. 

In addition, there is a mismatch of j profiles. The j distribution (or profile) 
within simulated halos has been found to scatter about a universal shape, with 
an excess of low-J (and high-j) material compared to the exponential disks ob- 
served [1, BD]. This mismatch has been demonstrated in an observed sample of 
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14 dwarf galaxies [18, BBS], which serves as the target for our modeling effort. 
BBS used for each halo the measured rotation curve and an assumed NFW pro- 
file to determine the halo virial quantities, with an average (V v i T ) ~ 60kms _1 . 
They then determined the baryonic spin parameter, averaging to (A^) <~ 0.07, 
significantly larger than the (A^ m ) ~ 0.035 of simulated halos, and then demon- 
strated the j-profile mismatch case by case. BBS also estimated the ratio of disk 
to DM mass to be (/a) ~ 0.04, about a factor of 3 smaller than the universal 
fraction, indicating significant gas loss. 

The spin catastrophe is commonly being associated with "over-cooling" , that 
the gas rapidly cools and becomes tightly bound in small halos. When such a 
halo spirals into a bigger halo, the baryonic component survives intact all the 
way to the center and thus transfers all its orbital j to the DM. It has therefore 
been speculated that energy feedback from supernova may remedy the problem 
by balancing the early cooling [e.g. 9]. A key idea is that the spin segregation 
between baryons and DM can go either way. While gas cooling tends to lower the 
baryonic spin, heating due to feedback would reduce this effect, and gas removal 
from small halos would even lead to higher baryonic spin. However, a realistic 
implementation of feedback has proved challenging [e.g. 17]. The feedback pro- 
cess has not yet been studied or implemented in satisfactory detail. We do not 
know yet whether they can indeed solve the CDM problems, and how. This moti- 
vates our attempt to first understand how the feedback scenario may work using 
a very simple semi-analytic model. Knowing that in a hierarchical scenario the 
halo fromation can be largely interpreted as a sequence of mergers, our model is 
based on a simple algorithm for the buildup of halo spin by adding up the orbital 
angular momenta of merging satellites [13, MDS; 19]. It matches well the spin 
distribution among halos in N-body simulations as well as the j profile within 
halos. This makes it a useful tool for understanding the over-cooling origin of 
the spin problem and for the attempt to cure it via feedback effects. Our work 
is described in more detail in [12, MD]. 

2 Buildup of Halo Spin by Mergers 

We characterize the angular momentum J of a galaxy by the modified spin 
parameter bull:01 A' = (J/M)/(\/2V v i I R v i T ). This quantity, which equals the 
standard A for an isothermal sphere and for an NFW profile with concentration 
c ~ 5, is straightforward to compute separately for the baryons and for the 
DM, X' h and A^ m . The distribution of In A' in the simulations is normal, with an 
average corresponding to Aq ~ 0.035 (compared to Ao — 0.042) and a standard 
deviation ay ~ 0.5. The "orbital- merger" model of MDS reproduces this spin 
distribution. To materialize this model we generate many random realizations 
of merger histories based on the Extended Press Schechter formalism [11] with 
slight adjustments, and for each merger tree we create random realizations of the 
orbital J added in each merger. The encounter parameters arc taken to mimic 
typical mergers, with the directions of the orbits drawn at random (or fine-tuned 
for a slight correlation between successive mergers as seen in simulations [5,16]. 
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The resultant distribution of halo spins matches the log-normal distribution 
obtained in the simulations. 

The cumulative mass distribution of j within simulated halos is fit by the 
universal function M(< j) = M v /j,j/(jo + j), with yu > 1 and j < j max = 
jo/(/i — 1) (BD). This is a simple power law, M(< j) oc j, for at least half the 
mass, with a possible bend characterized by /i. The other parameter, jo or j max , 
can be replaced by A'. The distribution of /x is Gaussian in In (it — 1), with a 
mean —0.6. The model also recovers these simulated j profiles. We create an 
M(< j) profile for each of the EPS model realizations by dividing the mass 
growth of the halo into bins and assigning to each the corresponding orbital J. 
A sample of profiles produced by this procedure and the distribution of [i values 
are shown in Figs. 1 and 2 of MD, demonstrating the match with the simulation 
results of BD. The model also reproduces the insensitivity to halo mass and 
rcdshift. 

The successes of the simple model in recovering both the distribution of spins 
and the j profiles makes it a useful tool for studying the j buildup. A new feature 
revealed by the model, which provides an interesting clue, is that the final halo 
spin is predominantly determined by the last major merger, while the many 
smaller satellites come in at different directions and therefore tend to sum up to 
a low j. If small satellites would lose gas before they merge into the halo, then 
much of the galactic gas would originate in big satellites, the final gas would lack 
the low-j component, and the baryonic spin would end up higher than the DM. 

3 Reproducing the Baryonic Spin Loss 

We can understand the j loss of baryons via a simple model including gas cooling, 
dynamical friction and tidal stripping for how the orbital j is converted into halo 
spin. First, the dynamical friction exerted by the halo on the satellite brings the 
satellite towards the halo center and thus transfers j from the orbit to the halo. 
Second, once satellite particles are tidally stripped they retain their j at the 
stripping point and add it directly to the halo. 

The mass loss at halo radius r can be estimated by evaluating the tidal 
radius £ t of the satellite at r via the resonance condition, m(£ t )/£t — M(r)/r 3 , 
where m(£) and M(r) are the mass profiles of the satellite and halo. If these two 
are self-similar, then this implies = r/i? v ; r , and the bound mass of the 

satellite is m[£ t (r)] oc M(r). A more accurate recipe for tidal stripping, tested 
with merger simulations, reveals that this is a good approximation in general 
[6], so we adopt it in our model. 

When exploring the effect of cooling, one assumes that initially the baryons 
follow the DM. As the gas cools, it contracts to a more compact configuration of 
radius i?b < Rdm- This spatial segregation in the satellite implies that the j-rich 
mass stripped at the early stages of the merger in the outer halo is dominated 
by DM, while the more compact baryons penetrate into the inner halo and lose 
more of their j via dynamical friction. The result is a net spin transfer from the 
baryons to the DM. Using the stripping recipe m(r) oc M(r), we obtain for the 
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final baryonic spin Jb/Jdm = {Rb/Rdm)- In the case of maximum cooling, the 
baryons dominate the halo center, = fbRdm, where /b — 0.13 is the universal 
baryon fraction. Fig. |l] shows the resultant baryonic spin distribution according 
to this model; there is a shift down to X' Q = 0.005, reproducing the spin crisis. 
The role of feedback would be to delay the cooling, increase i?b, and thus reduce 
the baryonic spin loss. 



V 



40 






. Over-Cooling 


Feedback 


30 




V„ = 0.040 - 






C,=0.43 ; 


20 




/ rL X\ = 0.068 










Jr \\ c,=o.6i 


10 


















0.00 0.05 0.10 0.00 0.05 0.10 0.15 

X 

Fig. 1. The effects of over-cooling and feedback on the spin distribution of baryons 
compared to the DM, for dwarf and bright galaxies. Log-normal fits are shown, with 
the mean and scatter quoted. Left: \' h is shifted down by an order of magnitude 
compared to A^ m , reproducing the spin catastrophe. Right: X' h in bright galaxies is 
boosted up by heating and partial blowout in incoming halos and roughly matches A^ m , 
while in dwarf galaxies \' h is boosted up further by the blowout in small satellites. 



4 Feedback 

Our approach here is to avoid the details of star formation and feedback and 
rather appeal to a very simple prescription for the effect of feedback as a function 
of the satellite's virial velocity, V sa t- Following the analysis of Dekel & Silk [8], 
we assume that the feedback by supernova-driven winds pumps energy into the 
gas and heats it uniformly to a temperature corresponding to a characteristic 
velocity Vfb, on the order of 100 km s -1 . We therefore assume that the spatial 
extent of the baryons is determined by the ratio Vfb/V Bai f The limit V sa t 3> Vfb, 
of massive, deep potential wells, corresponds to maximum cooling, i?b *C i?dm- 
In smaller halos, Still with V^ a t — Vfb, we expect the heating to balance the 
cooling and yield Rb — Rdm- Our model is therefore an interpolation between 
these limits, i?b = (Vfb/14at) 7l ^dm, with 71 an arbitrary exponent, which we set 
for now to be unity. 

If Vfb is larger than V^ a t, the feedback can cause gas blowout. We assume 
that partial blowout starts occurring once Vfb exceeds the escape velocity of 
the satellite, ~ V2V sa t: while total blowout is expected for Vfb 3> V, a t- We 
therefore parameterize the amount of gas that remains in the halo by another 
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interpolation, /d = (Vfb/Ks a t/\/2) 72 , with 72 an arbitrary exponent tentatively 
set to unity. We report here the results for the simplest choice 71 = 72 = 1, and 
explore the robustness of the results to different choices of 71 and 72 in MD. 

In Fig. |l| we demonstrate the effects of this feedback scheme, with Vfb — 
95kms _1 , on the distribution of X' h . We do it for two kinds of final halos, with 
V vlr = 60 and 220 km s^ 1 , representing dwarf and bright galaxies. The baryons 
in the bright galaxies end up with spins comparable to their DM halos, with 
Aq = 0.042, while in dwarfs they have significantly higher spins, with Aq = 
0.067. We learn that A b in dwarfs, which are build up by small satellites, is 
dominated by the blowout from these satellites, and it ends up with X' h > A dm . 
For bigger galaxies, which are largely made of bigger satellites, the dominant 
effect is the heating, with some contribution from blowout, together leading to 
a A b distribution similar to A dm , in general agreement with observations. 

In MD we explore a range of values for the exponents 71 and 72. For each 
choice we determine Vfb such that for the dwarfs (/d) = 0.04 as in BBS. We find 
that our results for dwarfs remain practically unchanged, while the results for 
bright galaxies have a weak dependence on the value of 71 in the range (0.5, 3). 

5 Model versus Observations 
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Fig. 2. Model realizations (with Vfb = 95kms _1 ) versus BBS observations of dwarf 
galaxies (shaded). Distribution of baryon fraction /d, spin parameter A', and j'-profile 
shape parameter £. Left: The model prediction for V v - a = 60kms _1 dwarfs, with 
significant blowout, is in agreement with the BBS data, while the bright galaxies retain 
most of their baryons. Middle: The predicted X' b distribution is in agreement with 
the dwarf data. Shown for comparison is the simulation result for dark halos, which is 
similar to that of bright galaxies. Right: The predicted distribution of £ for the baryons 
in dwarfs (heavy histogram) is shifted upwards compared to the DM (smooth curve) 
and the bright galaxies (light histogram), like the data, but its width is overestimated. 

The distribution of /d for the dwarfs observed by BBS is displayed in Fig. 0, 
showing values significantly lower than the universal value of /b — 0.13 and 
thus consistent with baryonic blowout. Shown for comparison are the model 
predictions for dwarf and bright galaxies. We enforced a match of the means for 
the dwarfs at (/d) = 0.04 by choosing Vfb = 95kms _1 , but the scatters are also 
in agreement. For bright galaxies, fd is typically lower than the universal value 
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by less than 50%, reflecting the limited fraction of small progenitors who lost 
their gas. 

Next, we compare predicted and observed spin distributions for dwarfs. We 
convert each value of A as quoted by BBS to A', and show their distribution 
in Fig. ||. the observed spins are significantly higher than the A' values of halos 
in cosmological simulations, with an average of X' ~ 0.07 compared to 0.035. 
Then shown is our model prediction with Vfb = 95kms _1 for the baryonic spin 
distribution in dwarfs. The effect of blowout brings the baryonic spins to a good 
agreement with the observed dwarfs. 

Fig. H shows the average j profiles and the scatter about them for the observed 
BBS dwarfs and for the corresponding model realizations compared to the typical 
j profile in halos by BD. We construct the baryonic j profile in each of our 
model realizations following the same method used to produce DM j profiles 
but now including feedback effects. The BBS dwarfs show low baryonic fractions 
(indicated by the integral under the histogram) and significant deficits of j at the 
two ends of the distribution compared to the halos. The profile for model dwarf 
galaxies is similar to the observations except for the very lowest j bin which is 
a 2cr overestimate, representing a spike in some of the model realizations. This 
spike may correspond to a low-j baryonic component that BBS fail to observe 
(faint halo stars?), or the spiky objects do not become disk dwarf galaxies, or our 
model needs to be improved. The high-j tail tends to be reduced in the baryons 
because it is often the result of a small satellite that comes in with its orbital J 
aligned with the halo spin, and now has lost its gas. 
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Fig. 3. Average j profiles (histogram) and the la scatter (shaded) for the observed 
dwarfs in comparison with the model dwarf and bright galaxies (Vfb = 95kms _1 ). The 
integral under the histogram is /d. Shown in comparison is the typical profile for DM 
halos by BD. 

Fig. H also shows the average model prediction for bright galaxies, Kir = 
220 km s -1 . They retain most of their baryons, so their profiles are less affected 
by blowout. The average model profile is in better agreement with an exponential 
disk, towards solving the j-profile discrepancy pointed our by BD. 
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A quantity used by BBS to characterize the shape of the j profile, as an 
alternative to the BD parameter fi, is £ = jtot/jmax- In Fig- § we also plot the 
distribution of this quantity in the BBS dataset in comparison with our model 
predictions for the baryons in dwarf and bright galaxies. The predicted ( dis- 
tribution for dwarfs is shifted upwards compared to the halos and the bright 
galaxies, in qualitative agreement with the BBS data, but the width is overesti- 
mated. 

6 Conclusion 

We devised a simple model to address the j problems of galaxy formation within 
CDM. By adding up the orbital J in random realizations of merger histories, the 
model successfully reproduces the simulated distribution of spins among halos 
(MDS) and the distribution of j within halos (MD). A simple analysis of how the 
merger orbital j turns into a spin profile provides a clue for how feedback effects 
in the satellite can resolve the spin problems. The idea is that the effective size 
of the gas component within the incoming halo determines its tidal stripping 
position in the big halo and thus its final remaining baryonic spin after the 
merger. The finding that the low-j material originates in many minor mergers, 
that tend to cancel each other's J, provides the clue for a possible solution to the 
j-profile mismatch problem. The blowout of gas from small incoming halos, which 
is more pronounced in satellites of dwarfs, would eliminate the low-j baryons in 
the merger product and increase the spin parameter, as observed. 

The feedback effects, including heating and blowout, are modeled as a func- 
tion of halo virial velocity, with one free parameter — the characteristic velocity 
Vh> corresponding to the feedback energy from supernovae. To match the low 
baryonic fraction observed in dwarfs it has to be Vfb ~ 100 km s -1 , consistent 
with the theoretical predictions [8]. This leads to an agreement between the 
model predictions and the observed disks, for the distribution of baryonic spin 
among galaxies and the baryonic j distribution within galaxies, both dwarfs and 
bright galaxies. 

We attempt to resolve the problems within the successful cosmological frame- 
work of CDM, by appealing to inevitable feedback effects. Another approach is 
to appeal to the Warm Dark Matter (WDM) scenario, despite the fact that it 
requires fine-tuning of the particle mass to ~ 1 keV. The main feature of WDM 
is the suppression of small halos and the corresponding mergers. While an N- 
body simulation of WDM [2] indicates the same j properties of halos (the same 
properties can also be obtained as a general result of tidal-torque theory, see 
MDS), one expects the cooling to be less efficient in the absence of small halos, 
and thus the baryonic spin to be higher. However, the j profile is still expected 
to be a problem, and the weaker feedback effects in the absence of small halos 
may not be enough for resolving it. These issues are yet to be studied in hydro 
simulations of WDM. 

Feedback effects may also provide the cure to the missing dwarf problem in 
CDM, where the predicted number of dwarf halos is much larger than the ob- 
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served number of dwarf galaxies [3] . While the number of dwarfs is automatically 
suppressed in WDM, it seems that the inclusion of the minimum inevitable feed- 
back effects would reduce the predicted number of dwarfs to significantly below 
the observed number and thus be an overkill (J. Bullock, private comm.). Finally, 
we find [7] that the key elements of our toy model — the tidal effects in mergers 
and the feedback in small halos — are also very relevant in understanding and 
resolving the third problem of CDM, where the halos in simulations typically 
show steep cusps in their inner profiles [14], while observations indicate flat cores 
at least in some galaxies [4] . An analysis of tidal effects explains the inevitable 
formation of an asymptotic cusp as long as satellites continue penetrating into 
the halo center. Feedback effects may puff up small satellites, make them disrupt 
in the outer halo and thus allow a stable core. 

The success of our toy model in matching several independent observations 
indicates that it indeed captures the relevant elements of the complex processes 
involved, and in particular that feedback effects may indeed provide the cure to 
all three problems of galaxy formation in CDM. The next natural step should 
be to incorporate a more sophisticated feedback recipe into the model using 
semi-analytic models and then full-scale cosmological simulations. 

This research has been supported by the Israel Science Foundation grant 
546/98, by the US-Israel Binational Science Foundation grant 98-00217, and by 
the German-Israeli Science Foundation grant 1-629-62.14/1999. 



References 

1. Bullock J. S., Dekel A., Kolatt T. S., Kravtsov A. V., Klypin A. A., Porciani C. & 
Primack J. R. 2001, ApJ, 555, 240 (BD) 

2. Bullock J. S., Kravtsov A. V. & Colin P. 2001, ApJL, in press 

3. Bullock J. S., Kravtsov A. V. & Weinberg D. H. 2000, ApJ, 539, 517 

4. de Blok, W.J.G., McGaugh, S.S., Bosnia, A. & Rubin, V.C. 2001, ApJL, 552, L23 

5. Dekel, A., Bullock, J.S., Porciani, C, Kravtsov, A.V., Kolatt, T.S., Klypin, A.A., 
& Primack, J.R. 2001, in Galaxy Disks and Disk Galaxie s, eds. J.G. Funes S.J. & 



E.M. Corsini (ASP Conference Series) ( |astro-pli/0011002[ ) 

6. Dekel, A. & Devor, J. 2002, in preparation 

7. Dekel, A., Arad, I., et al. 2002, in preparation 

8. Dekel A. & Silk J., 1986 ApJ, 303, 39 

9. Efstatiou, G. 2000, MNRAS, 317 697 

10. Fall S. M. & Efstathiou G. 1980, MNRAS, 193, 189 

11. Lacey C. & Cole S. 1993, MNRAS, 262, 627 

12. Mailer A. H. & Dekel A. 2002, MNRAS, submitted (MD) 

13. Mailer A. H., Dekel A. & Somerville R. S., 2001, MNRAS, in press (|astro 



ph/0105168j ) (MDS) 

14. Navarro J. F., Frenk C. S. & White S. D. M. 1995, MNRAS, 275, 56 

15. Navarro J. F. & Steinmetz M., 2000, ApJ, 538, 477 



16. Porciani, C, Dekel, A., & Hoffman, Y. 2001, MNRAS, in press ( bstro-ph/0105165| ) 

17. Thacker R. J. & Couchman H. M. P. 2001, ApJL, 555, L17 

18. van den Bosch F. C, Burkert A. & Swaters R. A. 2001, MNRAS, 326, 1205 (BBS) 

19. Vitvitska M., Klypin A., Kravtsov A. V., Bullock J. S., Primack J. R. & Wechsler 
R. H. 2001, ApJ, in press 



